%--------------------------------------------------------------------------
%
% Ephemeris computation using variable-order Radau IIA integrator with
% step-size control
%
% Last modified:   2018/02/11   M. Mahooti
%
%--------------------------------------------------------------------------
function Eph = Ephemeris(Y0, N_Step, Step,t0)

options = rdpset('RelTol',1e-13,'AbsTol',1e-16);
[t,yout] = radau(@Accel,(t0:Step:t0+N_Step*Step),Y0,options);
Eph(:,1) = t;
Eph(:,2:7) = yout;

